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The /ci-moment of a quark's Sivers function is known to be related to the corresponding twist- 
three quark-gluon correlation function T qi F(x,x). The two functions have been extracted from 
data for single-spin asymmetries in semi-inclusive deep inelastic scattering and in single-inclusive 
^ . hadron production in pp collisions, respectively. Performing a consistent comparison of the extracted 

' functions, we find that they show a "sign mismatch" : while the magnitude of the functions is roughly 

consistent, the fcx-moment of the Sivers function has opposite sign from that of T q ,F(x, x), both for 
up and for down quarks. Barring any inconsistencies in our theoretical understanding of the Sivers 
functions and their process dependence, the implication of this mismatch is that either, the Sivers 
effect is not dominantly responsible for the observed single-spin asymmetries in pp collisions or, 
the current semi-inclusive lepton scattering data do not sufficiently constrain the fcx-moment of the 
quark Sivers functions. Both possibilities strengthen the case for further experimental investigations 
of single-spin asymmetries in high-energy pp and ep scattering. 
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I. INTRODUCTION 



Since the observation of surprisingly large single transverse spin asymmetries (SSAs) in p^p — > irX at Femilab in the 
1980s [lj, the exploration of the physics behind the observed SSAs has become a very active research branch in hadron 
\ physics, and has played an important role in our efforts to understand QCD and nucleon structure [1, H|. Defined as 
• An = (ct(s_i_) — c(— s ±)) I (< J ( S ±) + a (— s ±))i the ratio of the difference and the sum of the cross sections when the 
| hadron's spin vector s± is flipped, significant SSAs have by now been consistently observed in various experiments 
, at different collision energies. These include semi-inclusive hadron production at low transverse momentum Phi in 
^ ; deep-inelastic scattering, IN^ -> i'hX, by the HERMES Collaboration at DESY [4], COMPASS at CERN @, and 
CLAS at Jefferson Lab @, as well as inclusive single- hadron production at high Phi in hadron-hadron collisions, 
p^p -> hX, by the STAR, PHENIX, and BRAHMS collaborations at RHIC Q. The observed large size of SSAs in 
hadronic scattering initially presented a challenge for QCD theorists [8J. Later two complementary mechanisms were 
proposed to describe the measured SSAs, and both of them have been quite successful phenomenologically [§-[l7j]- 

One mechanism relies on the so-called transverse momentum dependent (TMD) factorization [18l - l25l ]. and describes 
the SSAs in terms of the spin-dependent part of TMD parton distribution functions (PDFs), known as the Sivers 
functions [26], or TMD fragmentation functions (FFs), known as the Collins functions [27| ■ This TMD factorization 
approach is suitable for evaluating the SSAs of scattering processes with two very different momentum scales, Q\ ^> 
Q2 ^ Aqczj. The larger scale Q\ is necessary for using perturbative QCD, while the lower scale Q2 makes the 
observable sensitive to the parton's transverse motion. For example, the SSAs of hadron production at low Phi in 
lepton-hadron semi-inclusive deep inelastic scattering (SIDIS) have this characteristic property: Q ^> Phi. ~ Aqcd, 
and can be studied within the TMD factorization approach. 

The other mechanism generalizes the successful leading-power QCD collinear factorization formalism to the next- 
to-leading power in the expansion in 1/Q, where Q is the large momentum transfer of the collision, and describes 
the SSAs in terms of twist-3 transverse-spin-dependent three-parton correlation functions [28W3H , or a combination 
of the transversity distribution and three-parton fragmentation functions [30M32l |. This so-called twist-3 collinear 
factorization approach is more relevant to the SSAs for processes in which all observed momentum transfers Q are 
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much larger than Aqcd- This applies, for example, to the SSAs of inclusive single hadron production at high Ph± in 
p^p collisions. Although the two mechanisms describe the SSAs in two very different kinematic domains, they were 
shown to be equivalent in the overlap region where they both apply, and they thus provide a unified QCD description 
for the SSAs 

One of the potentially important contributions to the SSAs is the Sivers effect, which is generated by the initial- and 
final-state interactions between the struck parton and the spectators or the remnant of the polarized hadron [2 11 ] . The 
interactions provide the necessary phase that leads to the non-vanishing SSAs. In the TMD factorization approach, 
the role of these interactions is accounted for by including the appropriate color gau ge l inks into the definition of 
the TMD parton distributions, whose spin-dependent part defines the Sivers functions [23l. l34l. [35j . Since the details 
of the initial- and final-state interactions depend on the color flow of the scattering process, the form of the gauge 
links including the phase of the interactions is process dependent. Since the gauge links are included in the definition 



of the TMD parton distributions, the Sivers functions, too, are found to be process dependent [24|. Due to parity 
and time-reversal invariance of the strong interactions, the process dependence of the Sivers functions is effectively 
reduced to a sign change between their definitions in SIDIS and in Drell-Yan lepton-pair production in p^p collisions 
[23L f35j | . The predictive power of the TMD factorization approach relies on this modified universality of the Sivers 
functions. 

On the other hand, in the twist-3 collinear factorization approach, the process dependence of the initial- and final- 
state interactions is absorbed into the short-distance perturbative hard-part functions, while keeping the relevant twist- 
3 three-parton correlation functions universal or process- independent. The necessary phase for generating the non- 
vanishing SSAs arises from the quantum interference between a scattering amplitude with one active collinear parton 
and an amplitude with two active collinear partons. The SSAs are therefore proportional to the non-probabilistic 
three-parton correlation functions. Unlike the TMD parton distributions, which at given transverse momentum 
provide direct information on a parton's transverse motion, the twist-3 three-parton correlation functions provide a 
net asymmetry of the parton's transverse motion, after integration over all values of the parton's transverse momentum. 
As a result, the twist-3 three-parton correlation functions have a close connection with the transverse momentum k±- 
momcnt of TMD parton distributions. More precisely, the twist-3 quark-gluon correlation function, T qi p{x, x), often 
referred to as Efremov-Teryaev-Qiu-Sterman (ETQS) function, is equal to the first fc^-moment of the quark Sivers 
function f^(x, kj_) probed in SIDIS (or Drell-Yan) processes (23l. l33l. I3H] . 

Following the tremendous progress in experimental measurements of SSAs in recent years, the quark Sivers functions 
and the ETQS functions have been extracted for various quark flavors from the single-spin asymmetries in SIDIS and 
in pp scattering, respectively. In this paper, we examine the existing parameterizations of these two functions to see 
whether the first /c^-moments of Sivers functions are consistent with the existing twist-3 ETQS functions. Taking 
the quark Sivers functions f^p(x, kj_) extracted from SIDIS [ToL [Tlj. we evaluate their first fcj_-moments, and derive 
the ETQS functions T q ^p(x, x) with the help of the operator relation between the two functions [23, 33, 36]. We then 
compare the resulting "indirectly" obtained quark-gluon correlation functions with those "directly" extracted from 
the global fit [Tij to the SSAs for inclusive single hadron production in p^p collisions. In d oing so, we first observe that 
the sign convention adopted for the SSA in p^p — > hX in the previous literature [3, HH, l30j ] is in fact not consistent 
with that used for the experimental data. As a result, the signs of the T g> p(x, x) functions extracted in [14| need to be 
reversed. After this adjustment, we find that the twist-3 correlation functions T qt p(x,x) obtained in the two different 
ways have conflicting signs. 

The rest of our paper is organized as follows. In the next section, we briefly review the definitions of the quark 
Sivers functions and the twist-3 quark-gluon correlation functions (or ETQS functions). We recall the operator relation 
between the fcj_-moment of the Sivers functions and the ETQS functions, and discuss its limitations and the corrections 
to it. In Sec. IIII1 we present our findings regarding the sign "mismatch" between the existing parameterizations of 
quark Sivers functions and twist-3 quark-gluon correlation functions. We discuss the possible origins of this mismatch, 
and potential remedies. We also address the implications for phenomenology and propose further measurements to 
test the mechanisms for SSAs in hadronic processes. Finally, we give our conclusions and summary in Sec. IIVI An 
Appendix describes the derivation of the correct signs of the ETQS functions in single-inclusive hadron production in 
pp scattering. 



II. THE SIVERS FUNCTIONS AND THE ETQS FUNCTIONS 

In this section, we recall the definitions and relations between the quark Sivers functions and the twist-3 quark-gluon 
correlation functions, or ETQS functions. We use light-cone coordinates with the two light-like vectors 



n" = [1+0- Ox], »" = [0+1- Ox], 
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to project out the light-cone components: 



v 1 * g^u n- 



and v 



Qthvn v of any four- vector v 1 *. For the fully 



antisymmetric tensor e A " yp<T , we adopt the convention e 0123 = 1. We choose a frame in which the momentum of the 
transversely polarized hadron, p, is in the direction, with no transverse components: p p 



p + n fM . 



The quark Sivers functions for SIDIS kinematics with the transversely polarized proton moving in the +z-direction 
is defined through the following quark-field correlator [37j 
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is the gauge link consistent with the SIDIS process and V indicates path ordering [3J]. We note that the gauge link 
depends on the transverse separation £ i of the two field operators, which is responsible for the process dependence 
of TMD parton distribution functions J35|. Here it is worth pointing out that a different sign convention for the 
strong coupling constant g (for the interaction between the quark and the gluon) would lead to a different sign in the 
exponent of the gauge link in Eq. ^ (i.e., from ig to — ig). For deriving Eq. © we adopted the convention of the 



covariant derivative D a as 



(4) 



for the relevant part of the QCD Lagrangian density L — 'ipij^D^ip. Different conventions for (such as = 
— igAp) exist in the literature and in textbooks. Different conventions usually do not introduce any difference 
for a cross section, which has an even power in g. However, the single transverse spin asymmetry is proportional 
to the difference of two cross sections with the spin flipped, and the asymmetry is a consequence of an interference 
between scattering amplitudes of different phases, which is linearly proportional to ig. Therefore, one has to use the 
convention consistently in the theoretical definition and calculation of the single transverse spin asymmetry. 
Following the so-called Trento convention [38|], the correlator A4(x, k±) can be expanded as 



M(x, fcj 



/i(a:,fcj_)$- 
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where M is the nucleon mass, fi(x, k\) is the spin-averag ed TMD PDF, and ff T {x, k\) is the quark Sivers function. 
We note that a different convention for the Sivers function is commonly adopted in the phenomenological studies by 
the Torino group PUSH- Here a function A N f q / A t [x, k±) is introduced which is defined from 



f g /At(x,kx) = Tr 



~ij,M(x,kx) 



1 



= h(x,ki) + -A»f q/At ( 



,kx)sx ■ (p x kj.). 



(6) 



The relation between A f q /At an d t ne Sivers function in the Trento convention is 

2fc ± 



A JV / g /At(z,fc±) 



M 



/it ( x > k*x)- 



In the twist-3 collinear factorization approach, the ETQS function T q ,p(x,x) is defined as [lil [l6| 



T q<F (x,x) 



4n 



where Vjo,f] and Vjf^j are the gauge links along the "— " light-cone direction and are given by 



ig l dr/ A + (i] ) 



(7) 



(8) 



(9) 



Here the sign convention for coupling constant g is the same as that in Eq. (j4|). Choice of, for example, the convention 
with Dp = dp — igAp would change the sign of the exponent. Within the collinear factorization approach, it is assumed 
that the typical transverse momentum of all active partons is much smaller than the hard scale of the scattering 
process, Q. Up to power corrections in 1/Q, the transverse momenta of all active partons are completely integrated 
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into non-perturbative PDFs, FFs, or the correlation functions. Consequently, unlike for the TMD distributions, all 
field operators defining the non-perturbative functions in the collinear factorization approach are evaluated at the 
same light-cone separation with zero "+" and "_L" components, as shown for example in Eq. (JHJ. 

Since the quark-gluon correlation functions in the collinear factorization approach have all their active partons' 
transverse momenta integrated, these correlation functions can be related to fcj_-moments of the TMD parton distri- 
bution functions. It was shown at the operator level [23|, [U [3(| that the ETQS function T^f{x,x) is closely related 
to the fci-moment of Sivers function: 

gT q;F (x,x) = - J ftxM/i'f^Jm (10) 

where the subscript "SIDIS" emphasizes that the Sivers functions here are probed in the SIDIS process. We stress 
again the importance of the sign convention for the coupling constant g in the definition of the gauge link. If the sign 
convention used to define T qi F(x, x) is different from that in the definition of f^p(x, k\), the difference will introduce 
an extra factor "—1" in the relation between these two functions, so that there will be no minus sign on the right-hand 
side of Eq. (pTT|) . 

We emphasize that the operator definition in Eq. ([5]) does not completely fix the quark-gluon correlation function 
T qi p(x, x), unless the renormalization scheme is specified. As is well known from the case of ordinary PDFs, the matrix 
element in Eq. (JSJ) is ultraviolet (UV) divergent [Hj]. Like in the case of PDFs, the quark-gluon correlation function 
is really defined in terms of the QCD factorization formalism. The leading UV divergent (the large k±) region of the 
matrix element on the right- hand- side of Eq. corresponds to the region of phase space with large parton virtuality, 
and is required by factorization to be moved from the matrix element into the perturbatively calculated short-distance 
functions. The removal or subtraction of the UV divergence is not unique, which leads to the factorization scheme 
and scale (fj,) dependence of the correlation functions T q ^p(x, x, /i) 40]. In this way, also the relation in Eq. (|10l) is 
subject to the UV subtractions and the adopted factorization scheme, and hence not a unique identity. That said, the 
relation (TTU1) provides a natural "zeroth-order" connection between the Sivers and the ETQS functions. It plays an 
important role in establishing the consistency between the TMD factorization approach and the collinear twist-three 
quark-gluon correlation approach in the descriptions of the SSAs in SIDIS and the Drell-Yan process [33(. It also is a 
useful starting point for phenomenological studies and is of much help in testing the various constraints on the quark 
Sivers and quark-gluon correlation functions. In the following, we will therefore make use of relation (1101) . keeping 
however in mind the caveats we have made regarding UV renormalization. 



III. THE "SIGN MISMATCH" 

The quark Sivers functions f^(x,k^_) (or equivalently, A N f q / A t (x, k±)) and the twist-3 quark-gluon correlation 
functions T q} p(x,x) have been extracted from experimental data on SSAs for single hadron production in SIDIS and 
in hadron-hadron scattering, respectively. In this section, we compare the existing parameterizations of these two 
functions and present our findings concerning the "sign mismatch" . We also introduce and discuss various loopholes 
that might resolve the apparent inconsistency. 

So far the quark Sivers functions have been extracted from the A s ^^ h azimuthal asymmetries in SIDIS. We 
consider two such parametrizations here. One is from Ref. [k} (we refer it as "old Sivers"), the other one ("new Sivers") 
from Ref. [ll| • They both parametrize the spin-averaged TMD PDFs ff (x, kj_) and Sivers functions A N f q / h t (x, k±) 
for each quark flavor q in the form 

f?(x,k 2 ± ) = mx)g(k ± ), (11) 
A N f q/hf (x,k ± ) = 2Af q (x)f?(x)h(k ± )g(k ± ), (12) 

where ff (x) is the quark's spin-averaged collinear PDF, Af q (x) is a fitted function whose functional form is not relevant 
for our discussion below, and g(k±) is assumed to have a Gaussian form, 

= ^2T e ^ i/<fci> ( 13 ) 

with a fitting parameter (k^_) for the width. However, the two parameterizations adopt different functional forms for 
the &j_-dependence of the Sivers function: 

old Sivers: h(k±) = , 2 ? fc±M ° , (14) 
k\ + Mq y ' 

new Sivers: h(k±) = yfYe^-e-^l™* , (15) 
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where Mq and Mi are fitted parameters. 

Since for both parameterizations the fc^-dependence is assumed to be decoupled from the x-dependence, we can 
derive the cc-dependence of the associated twist-3 quark-gluon correlation T q ^p{x, x) analytically, using the relation in 
Eq. (|10l) . By substituting the parameterization of the Sivers function in Eq. (fT2")) into the right- hand-side of Eq. (|10p , 
and using the fitting parameters extracted in Refs. (lOj and fllj ]. we obtain the following two parameterizations for 
the correlation function T qt p{x,x): 

gTq t p(x, x) | i d Sivors = 0.40/f(a;) J V g (a;)|oid I (16) 
gT qyF {x,x)\ 

new Sivcrs new • 

(17) 

From the existing data, the best constrained Sivers functions are those of u and d quarks. Using the fitted functions 
N q {x)\ \& and M q {x)\ ncw from Refs. [l(| and [llj], respectively, we plot the drived quark-gluon correlation functions 
x gT Ut p(x 7 x) (left) and x gTd t F(x,x) (right) in Fig. [I] The dashed lines are for the quark-gluon correlation functions 
obtained by using the new Sivers parameterization, while the dotted lines arc for the old Sivers parameterization. We 
find that for these "indirectly" obtained quark-gluon correlation functions, T u ^f{x,x) is positive, while Td,F{x, x) is 
negative. 




FIG. 1: The quark-gluon correlation function gT qt F(x, x) as a function of momentum fraction x for it-quarks (left) and d-quarks 
(right). The dashed (dotted) lines are gT q: F(x, x)| now sivcrs {gT q ^{x, x)\ \& sivcrs) obtained by taking the fc^-moments of the 
corresponding quark Sivers functions according to the right-hand-side of Eq. (|10p . The solid lines represent the correlation 
functions extracted directly from data on SSAs for inclusive pion production in proton-proton collisions, p^p — >■ n + X [l4| ]. 
after correcting for the sign convention (see text). 



On the other hand, the ETQS function T qy p(x,x) can be "directly" extracted from data on SSAs for inclusive 
single hadron production in hadronic collisions, p^p — > h(Ph±, y) + X, assuming these asymmetries are predominantly 
generated by the Sivers effect (or rather, its twist-3 counterpart). Such SSAs have been measured at sufficiently 
large transverse momentum P hl _ by the E704 Collaboration at Fermilab [J, and the STAR, PHENIX, and BRAHMS 
collaborations at RHIC 7] • Since they depend only on one large momentum scale Ph± , these SSAs are better studied 
in the collinear factorization approach, where they may be generated by three possible mechanisms: (1) the twist- 
3 quark-gluon and tri-gluon correlation functions of the polarized hadron, (2) the transversity distribution of the 
polarized hadron combined with the twist-3 quark-gluon fragmentation functions to the observed hadron, and (3) the 
transversity distribution combined with possible twist-3 unpolarized quark-gluon correlation functions j3(j. It was 



found that the third mechanism only makes a small contribution [31| . By assuming that the observed SSAs are mainly 



generated by the ETQS functions T qt p(x, x), a set of T qt p(x, x) was extracted by a global fitting procedure [14[. In the 
course of our investigations, we have revisited the sign convention adopted in [14| . which was based on the earlier work 
0, HH, We have discovered that the convention was at odds with that chosen in the experimental studies. The 
inconsistency can be traced to the value of the contracted Levi-Civita tensor appropriate for the spin asymmetry. We 
provide a detailed discussion of this issue in the Appendix. Correctin g th e sign convention of 0, H^, Ho] means that 
one needs to change the signs of the T q ^p(x,x) functions extracted in [lj]. We plot the resulting "directly" extracted 
T q ^p{x,x) functions as solid lines in Fig. [TJ along with the previous ones derived "indirectly" from the fc^-moment of 
the quark Sivers functions. Surprisingly, we find that the two sets of functions have opposite signs, both for up and 
for down quarks. 
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At first sight, it may seem that we have created a problem where none used to be. After all, the sign mismatch we 
find becomes apparent only after we have changed the signs of the T q< p (x, x) functions of [14]]. However, the basic 
problem is easy to see: as we discussed in the Introduction, the Sivers contributions to the single spin asymmetries 
depend on initial or final state interactions in the scattering processes. The SSA in SIDIS comes from a final state 
interaction. A negative up-quark Sivers function is known to generate a positive SIDIS spin asymmetry for n + 
production. In p^p — » tt + X at forward rapidities, however, the main partonic channel is ug — » ug, for which initial- 
state interactions play the dominant role, resulting in negative partonic hard-scattering functions. Therefore, if the 
Sivers mechanism (or its twist-3 variant via Eq. (|10p ) is primarily responsible for the SSA in this process, one would 
expect a negative asymmetry for 7r + , contrary to what is observed. Thus the T q ^p(x, x) functions needed to describe 
the RHIC single-spin asymmetries cannot have the signs suggested by Eq. (fTU|) . We note that in these considerations, 
one has to carefully take into account the experimental definitions of the SSAs; see the Appendix for some details. 

There are two main caveats regarding the sign mismatch. The first one is that the integral over k± in Eq. (fTOj) 
might produce a different sign from that of f^(x, k\) itself in the region of k± where it is constrained by data. The 
HERMES SIDIS data that are mostly relevant for the extraction of fyS(x, kj_) are at a relatively modest Q 2 ~ 2.4 
GeV 2 . Since the TMD factorization formalism is valid only for k± <C Q, the data constrain the function and its sign 
only at very low k± ~ Aqcd • The existing parameterizations of the quark Sivers functions [ToL ITlT ] assume a purely 
Gaussian form of the fcj_-dependence and hence would not allow a sign change of the function at some This 
leads to significant uncertainties in the determination of the twist-3 quark-gluon correlation functions via Eq. (jlOl) . 
because taking the k± moment enhances the contribution from the unknown larger-fc^ region. Also, the issue of UV 
renormalization discussed in the previous section becomes relevant for the k± moment. We note that future SIDIS 
experiments at an Electron Ion Collider would have the kinematic reach to precisely map out the k± dependence, 
and to allow measurements of the transverse-momentum weighted asymmetries, providing direct access to the twist-3 
quark-gluon correlation functions. In this way, reliable comparisons with the correlation functions extracted from pp 
collisions would become possible. 

It is worth keeping in mind that the SIDIS and pp single-spin asymmetry data also probe slightly different values 
of x. The former reach up to x ~ 0.4, while the latter mostly access yet larger values, x ~ 0.6. While it is in principle 
possible that a rapid sign change could occur towards large x which would explain the mismatch, there is nothing in 
the SIDIS data or the pp data with a sufficiently large xp coverage that would indicate such a behavior, and we do 
not consider this to be a likely scenario. 

The second possibility is that there are other significant contributions to the SSAs for single hadron production 
in p^p collisions, besides the Sivers mechanism. In addition to the asymmetry due to the spin-dependent twist-3 
quark-gluon correlation functions of the polarized hadron, the SSAs in hadronic collisions may also be generated at 
the hadronization stage by a combination of the transversity distribution of the polarized hadron and the twist-3 
quark-gluon fragmentation functions [32( |. which is effectively a representation of the Collins effect in the collinear 
factorization approach. If this mechanism makes a large contribution to the hadronic SSAs, with sign opposite to 
that by the T q ^p(x,x), it might explain the observed features. Unlike the measurement of SSAs in SIDIS, where 
the Sivers effect and the Collins effect can be separated by using different azimuthal angle weighting, the two effects 
cannot be separated in single-hadron inclusive production in hadronic collisions. Nonetheless, other measurements 
are available in pp scattering that would allow to disentangle them. The prime example is the Drell-Yan process, 
which allows direct access to the Sivers or T q ^p(x, x) functions [23l. [35|. A similar role could be played by photon pair 
production [44| . Here we will briefly consider two further processes that have the advantage of being somewhat more 
copious at RHIC. 

In order to get clean access to the quark-gluon correlation functions T q> p(x, x), we need observables that are 
not sensitive to the details of the hadronization stage. At RHIC, for example, direct photon production [25[ at 
large transverse momentum Ph± and inclusive single jet production at large transverse jet energy are two promising 
observables of this kind. Since the fragmentation contribution to prompt photon production at large px is much 
smaller than the direct contribution at RHIC energies, in particular if photon isolation cuts are imposed, the SSAs of 
these two observables could provide direct information on the twist-3 quark-gluon correlation functions, thus allowing 
to see if their signs are consistent with those derived from Eq. (| 10[) . 

In Fig. [2] we present our estimates for the SSAs for direct photon (left) and inclusive single jet production (right) 
in p^p collisions at \J~S = 200 GeV. We consider production in the forward region of the polarized proton, so that 
to a good approximation we only need to include the valen ce q uark contribution for the polarized beam. For the 
relevant unpolarized PDFs, we use those specified in Ref. fiol llll llij] correspondingly. The solid curves represent the 
SSAs calculated by using the directly extracted T q ^p(x, x) (with T Ut p(x, x) < and T^p{x, x) > 0), which were shown 
as the solid lines in Fig. [TJ The dashed and dotted curves show the SSAs calculated by using the indirectly derived 
T qi p(x,x) from Eq. (|16[) and Eq. (|17[) . respectively, which again were shown by the same line patterns in Fig. [T] and 
have T Uj p(x,x) > and T^p[x,x) < 0. The results in Fig. [5] demonstrate that positive An for direct photon and 
inclusive single jet production should be expected at RHIC if the "directly" extracted T q ,p(x,x) are correct. If, on 
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the other hand, the signs of T g> p(x,x) follow Eqs. (|T6|) and (flTl) . negative values for the An for the two processes are 
predicted. We note that direct photon and inclusive single jet production both receive contributions from the u and 
d quark ETQS functions. Since these have opposite signs and rather similar magnitude, their effects cancel to some 
degree for jet production. For photons, the situation is more favorable thanks to the weighting by the quark's charge 
squared, which explains why here the spin asymmetries are overall larger. 
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FIG. 2: The SSAs for direct photon (left) and single inclusive jet (right) production in p^p collisions at \/~S = 200 GeV, as 
functions of xf for rapidity y = 3.3. The various curves correspond to the T q< F(x,x) shown in Fig. [1] 



IV. SUMMARY 



We have computed the /c^-moments for two parameterizations of up and down quark Sivers functions determined 
from semi- inclusive lepton scattering data given in [Tol [Tl| . These are related to the quark-gluon correlation functions 
T q< p(x,x) relevant for the description of single-spin asymmetries in single hadron production in pp scattering. The 
latter have in the past been extracted from RHIC data Correcting an inconsistency in previous theoretical 

treatments of the spin asymmetries in pp scattering, we have found that the resulting T q p(x, x) functions have signs 
opposite to those predicted from the analysis of the fcj_-moments of the Sivers functions. We have discussed various 
possible explanations for this apparent discrepancy. 

Our finding highlights the importance of additional measurements of single-spin asymmetries. Measurements of the 
k± dependence of the Sivers functions with wide kinematic reach would be feasible at an Electron Ion Collider and 
should shed light on the contributions from various fcj_-regions to the moment of the Sivers functions. We have also 
shown that Am measurements for jet and direct photon production in pp collisions at RHIC should be valuable tools 
for a cleaner determination of the quark-gluon correlation functions T qt p(x,x). 
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Appendix: The sign of T q] F(x,x) in inclusive hadron production 

In this appendix, we demonstrate why the SSA data for p^p — > hX require T Uj p(x,x) < and Td t p(x, x) > 0, if 
the ETQS functions are the dominant sources of the observed asymmetries. 

We start with the QCD factorization formalism for the spin-averaged cross section for inclusive single particle 
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production in hadronic collisions, A^(S±) + B — >• h(Ph±) + X: 

Eh ^P h = f £ / / f^/aCO / f /.M(x)^ c (S,t,fi)«(S + * + fi), (A.l) 

a,6.c 

where f a /A{%) and /t/s(x') are the PDFs, D c ^h(z) are the FFs, and H u b ^ c are the partonic hard-scattering functions, 
with s, i, and u the Mandelstam variables at the parton level. Including only the contributions by the twist-3 quark- 
gluon correlation functions, the spin-dependent cross section dAa(s±) = [dcr(sj_) — da(— sj_)]/2 is given by 

dAcr(.s ± ) /- d2 , f r dx! r , f dx n / £ ^ Si " ft \ 

iJ ab ^ c (s,f,M)(5(s + £ + u) , (A.2) 

where the relevant hard-scattering functions H a b^ c (s,t,u) can be written as 

H ab ^ c (s, t, fi) - ff^ c (s, t, fi) + fl£_, c (a, t, fi) + j) . ( A - 3 ) 

with ff£ 6 _> c and H a F h ^ tc representing the contributions from initial- and final-state interactions, respectively. The 
explicit forms of H^ b ^ c , H^ b ^ c , and H^ b ^ c are given in (l4j. It is important to point out that the spin-dependent 
cross section in Eq. (|A.2[) is calculated from an interference between two partonic amplitudes. It thus depends on 
the sign convention for the coupling constant g; the form given in Eq. (|A.2[) is based on the convention in Eq. 
If one uses the other sign convention for the covariant derivative, there will be an extra minus sign appearing on the 
right-hand side of Eq. (|A.2|) . which would be compensated by an extra sign in Eq. (ITU1) . 
The SSA, A n, is given by the ratio of spin-dependent and spin-averaged cross sections: 

„ dAats i ) / „ da , , . , . s 

E h d3 p f ' / E h-#[p- h = sin(0 s - </> h ), (A.4) 

where 4>h and 4> s are the azimuthal angles of the hadron transverse momentum Phj_ and the spin vector s±, respectively. 
The absolute sign of An depends on the choice of frame and the coordinate system. In experiment the following 
convention is used: positive values of An correspond to a larger cross section for hadron production to the beam's left 
when the beam's proton spin is vertically upward .'!() , as sketched in Fig. IA.1I In the center-of-mass frame of A and 
B, a convenient coordinate system (consistent with the experimental convention) is given by choosing the polarized 
nucleon A to move along +z, the unpolarized B along —z, the spin vector sj_ along y, and the produced hadron's 
transverse momentum Phi. along the x-direction. In this frame, 4>h = 0; 0s = tt/2, and 

e Ph ^ nn = -\P h ±\\s ± \. (A.5) 

We note at this point that there is an overall sign error in [3(J and consequently in [31, because in these papers the 



T a F (x, x) - x—T aF (x, x) 
dx 




FIG. A.l: Illustration of the sign convention for An'- positive An means that more hadrons are produced to the left of the 
beam direction when the beam's spin is vertically upward. 



choice e Ph±s±nfi > was made (see Eq. (73) of [30], in contrast to Eq. (|A.5j) above). 

In the forward direction, qg — > qg is the dominant partonic scattering channel for inclusive single hadron production. 
The corresponding hard-scattering functions are given by 14 1 



H. 



2N? 



s 
u 



u 
s 



1 



2N^ su 
N? - 1 ¥ 



|t|<s~|u| 



"2s 2 



(A, 
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H, 



qg^qg 



2{N 2 - 1) 



2N 2 {N 2 _ 1) 



s u 




, U 2 ' 


\t | <C s ~ \u\ 






1 - N 2 ^- 
t 2 . 


u s 







2(N 2 1) 



2s 



21 



S 



1 + 2Nt 



, sit 



1 



N 2 -l 



t 2 

'2s 2 



t 2 



(A.7) 
(A.8) 



This shows that both H qg 

H gg ^ qg . Furthermore it is clear that the SSA in ir + production is mainly sensitive to T u ,p(x, x), while the one for it 



qg and Hqg^qg have opposite sign to that of the spin-averaged hard-scattering function 



production probes Td,F(x, x). Since 



>0, 



(A.9) 



we conclude from Eq. (|A.2[) that the observed positive SSAs for 7r + production indicates a negative T U; p{x,x), while 
the observed negative asymmetry for n~ production indicates a positive p(x,x), as shown by the solid curves in 
Fig.ffl 

To conclude this appendix, we demonstrate the apparent "sign mismatch" again numerically, by evaluating the 
SSAs for inclusive single hadron production using the ETQS functions indirectly derived via Eq. (1101) from the quark 
Sivers functions in Eqs. (fl"6|) and (fT7|) . The results are shown in Fig. IA.2I As expected, the signs of the calculated 
SSAs are opposite to those observed experimentally. 




FIG. A. 2: The SSA, An, for inclusive single pion production in p^p — > ty + X at yfs — 200 GeV, as a function of xf and at 
rapidity y = 3.7, evaluated by using the old Sivers functions in Eq. (fTo]) (left), and the new Sivers functions in Eq. (TITl) (right). 
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